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ABSTRACT 

We discuss the generalisation of the so-called traditional approximation, well 
known in geophysics, to general relativity. We show that the approximation is 
applicable for rotating relativistic stars provided that one focuses on relatively 
thin radial shells. This means that the framework can be used to study waves 
in neutron star oceans. We demonstrate that, once the effects of the relativistic 
frame-dragging are accounted for, the angular problem reduces to Laplace's 
tidal equation. We derive the dispersion relation for various classes of waves in 
a neutron star ocean and show that the combined effects of the frame-dragging 
and the gravitational redshift typically lower the frequency of a mode by about 
20%. 



1 MOTIVATION 



Accreting neutron stars exhibit a variety of quasiperiodic phenomena, with frequencies rang- 
ing from a few Hz to well above 1 kHz. Many of the observed oscillations are likely associated 
with processes in the accretion disk, but it is also pr obable that some of the observed fea- 



tures are linked to waves in the neutron star ocean (Bildsten and Cutler 



Strohmaver and Bildsten 



1995 



HeyJ 



2001 



2003). For example, it has been observed that neutron stars that 
accrete matter at rates > lO~ 9 M yr _1 (the so-called Z-sources) show quasi-periodic oscilla- 
tions with frequencies of ~ 6 Hz. It has also been argued that the Z-sources are covered by 
massive oceans comprising a mass of ~ 10~ 6 M^ composed of a degenerate liquid of light 
elements (C, O, Ne, Mg,..). This led lBildsten and Cutler! (|l995T ) to argue that the observed 



oscillations may correspond to ocean g-modes. 



To analyse this suggestion in more detail 



Bildsten. Ushomirskv and Cutler! (|l996l ) inves- 



tigated the nature of low-frequency waves in the ocean of a rotating neutron star. They 
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considered slowly rotating stars (far below the breakup limit) in Newtonian gravity, and 
made use of the assumption that Q ^> a, where Q is the rotation frequency and a is the 
mode frequency. This is a reasonable assumption since neutron stars accreting at a high 



rate a re ex pected to spin much faster than a few H z. The studies of 



Bildsten and Cutler 



(1995|) and Bildste n. Ushomirskv and Cutler! (|1996f ) both made use of the so-called "tra- 
ditional" approximation, which has its origins in geophysics. This approximation leads 
to a separation of variables for the linearised Euler equations, where the angular prob- 
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rors ( 


Miiller. Kellv and 0' Brien 


1993 


Lee and Saio 


1997: 


Townsend 


20C 


3). 



Mile; 



The purpose of this paper is to extend the "traditional approximat ion" to general rela 



Bildsten and Cutler 



tivity. The main motivation for this is tha t we know from th e studi es of 

1 1 irhi in 

(1995); Bil dsten. Ushomirskv and Cutlerl (1996) Townsend! 1)20031 ) that the approximation 
provides a relatively simple way of analysing ocean waves in rotating Newtonian stars. Yet 
there are several reasons for studying this problem in the framework of general relativ- 
ity. First of all, neutron stars are highly relativistic objects with GM/Rc 2 ~ 0.15 — 0.2. 
This means that one must account for the gravitational redshift when comparing observed 
oscillations to theoretical models. A second relativistic effect that needs to be accounted 
for is the rotational frame-dragging. Finally, the traditional approximation scheme is in- 
teresting from a conceptual point of view. Especially since it provides an alternative ap- 



proach to existing work on the low - 



reque ncy oscillations of rotating relativistic stars (see 



Lockitch. Friedman and Andersson! (J2003J) for a discussion and references to the relevant 
literature) . 



2 THE RELATIVISTIC EQUATIONS OF MOTION 

We begin by deriving the equations that govern a small amplitude ocean wave on a relativistic 
star. We consider a star which rotates uniformly at frequency £1 <C (GM/R 3 ) 1 ^ 2 , where 
M and R are the neutron star mass and radius. This means that the star rotates at a 
small fraction of the break-up velocity, allowing us to neglect the "centrifugal force" and 
treat the unperturbed star as spherical. This approximation is justified for most observed 
neutron stars. However, it is imperative to account for the Coriolis force. In particular since 
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it leads to the existence of a new class of oscillation modes, the so-called "inertial modes" 



( Provost. Berthomien and Rocca 



mi 



Lockitch. Friedman and Andersson 
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I 



The spacetime of a slowly rotating star is described by the metric 

ds 2 = -e v dt 2 + e x dr 2 - 2ur 2 sin 2 6dt dip + r 2 (d6 2 + sin 2 9d<p 2 ), (1) 

where u represents the dragging of inertial frames (|Hartldll967l ). Physically, uj represents the 
angular velocity of a "zero-angular momentum observer" (for which the associated specific 
angular momentum, u^, vanishes) relative to infinity. For an observer co-moving with the 
fluid, the four-velocity moves on wordlines of constant t and 9, and the corresponding angular 
velocity Q is defined as the ratio 



dt 



(2) 



Combining this with the requirement that 



u» Utl = -l (3) 

we have 

u» = [e-"/ 2 ,0,0, Qe~ u/2 }. (4) 

We wish to consider small perturbations of the star. Since we are mainly interested in 
ocean waves we will use the Cowling approximation; i.e. assume that all metric perturbations 
can be neglected. This is a reasonable assumption given the relatively low density of a neutron 
star ocean. One would, for example, not expect the ocean waves to radiate appreciable 
amounts of gravitational waves. We also need to neglect the dynamic nature of spacetime if 
we want to implement the "traditional approximation" in general relativity. 

In order to derive the equations that govern the fluid motions we assume that the stellar 
background is appropriately described by a perfect fluid equation of state with p = p(p), 
which means that the energy-momentum tensor takes the form 

= pu^u u + pq"", (5) 

where q^ v is the projection orthogonal to u^; 

q»» = + (6) 

and the scalars p and p denote the energy density and pressure, as measured by a co-moving 
observer (an observer with four- velocity respectively. 
The equations of motion follow from the conservation law 

V M T^ = 0, (7) 
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as well as the Einstein field equations. It is well-known that, in the case of slow-rotation, this 
problem reduces to the standard Tolman-Oppenheimer-Volkoff equations for a non-r otating 



star, plus a single ordinary differential equation determining the frame-dragging (jHartk 



19671 ) . These equations can be written 



y = 2e A (M-47rpr^ (g) 
, 2e A (M + 4vrpr 3 ) 

V = ~2 ( 9 ) 

/ P + P , / ln % 

where primes denote derivatives with respect to r, and M = M(r) is the "mass within radius 
r" , which follows from the relation 

=(!--) ■ (ID 
To incorporate the leading order rotational effects, we need to solve an additional equation 
that describes the frame-dragging, uo. It is convenient to use the difference Co = Q — u, corre- 
sponding to the rotation frequency as measured by local zero-angular momentum observers, 
which is determined by the equation 

[ e -(H-A)/2 + l e -(,+A)/2 -/ _ 167r(p + p)e (A-,)/2 - = g ? (12) 

We also need boundary conditions for uj at the origin and spatial infinity At the origin 
we require that Co is regular. In the vacuum outside the star, equation (fT^jl can be solved 
analytically, and we have 

2 7 

u = n-— , (13) 

where J is the total angular momentum of the star. This relation can be used to provide 
boundary conditions for Co (and its derivative) at the surface of the star in terms of Q and 
J. Specifically, the solution to equation (fT2^) is normalised by requiring that 

„ (R) + A^fl) = n (14) 

In order to derive the equations that govern perturbations, it is useful to work with 
projec tions of equation (J7J) along, and orthogonal to, the four- velocity, cf. llpser and Lindblom 
( 19921 ) . Projecting along m m we have the relativistic continuity equation 

^+(p + /))V/ = 0. (15) 

The orthogonal projection leads to the relativistic Euler equations 

(p + pKV^ + g£V M p = 0. (16) 
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Perturbing equation (|15p. assuming that the Cowling approximation is made (i.e. taking 
Sg^u = 0), we get 



Su^V^p + w M V M <5p + (p + p) V^u* = 0. 

Meanwhile, the perturbed Euler equations can be written as 

(5p + 5pKVX + (p + p)Su"V^u K + {g VK + V,<5p + (p + p)^V M 5w K 
- fcW(p + p)V^ = 0. 

The perturbed fluid velocity has (quite generally) components 

5m* = (n-uj)r 2 sin 2 de- v/2 H 2 e i{(Tt+mv \ 

5m 9 = e "/2jJ ie iM+w) ) 

(J^v = e " /2 fl2e <(a * +m *' ) , 



(17) 



(18) 



(19) 
(20) 
(21) 
(22) 



where <kt* has been fixed by the requirement that the overall four-velocity + <5w M is 
normalised in the standard way. Here W, Hi and Hi are functions of r and 9, and a is the 
oscillation frequency observed in the inertial frame. In addition we let 



Sp = Pl e i{at+mv) . 



(23) 
(24) 



where pi and pi are functions of r and 9. Given these definitions, the continuity equation 
(fT7l) leads to 



i(a + mQ)pi + (p + p)e u W + 

OHi cos 9 



e u W 



+ (p + pK 



09 sin 6* 



+ im(p + p)e v if 2 = 



(25) 



The r- component of (|18|) gives 



2 r (p + p) [r (w - to) v' - Jr + 2 (O - w)] sin 2 fl# 2 
-2 pi - (pi + pi) i/ - 2z((7 + mil) (p + p) e x W = 0, 

while the 6* and y> components lead to, respectively, 

2 (p + p)r 2 (fi-cj)cos#sin#iJ 2 - % -i(«7 + mfi) (p + p)r 2 #! = 0, 



(26) 



(27) 
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and 

2(p + p)r 2 (VL — uj) sm 9 cos 9 Hi + imp 1 + i(a + mVL)(p + p)r 2 sin 2 9 H 2 

+ {p + p)[(fl-u)(2r-r 2 u') - r 2 uj'}sm 2 9W = 0. (28) 

It is worth emphasising the following point: In deriving the above equations we have 
neglected terms proportional to afl in the (^-component of the continuity equation (fTHjl . Such 
an approximation is not required in the Newtonian the terms that we have omitted 

are post-Newtonian contributions. Specifically, the complete ^-component of equation ()18|) 
includes the term 

,cr(n-u) 2 . 2 
% - r sm 9p\ 

cr 

where we have reinstated the speed of light (c). This term vanishes in the Newtonian limit 
(c — > oo), but in the relativistic case we can neglect it only when a is suitably small. 

We do not, however, expect this difference between the Newtonian and the relativistic 
equations to have great effect on the ocean waves that we are interested in. The main reason 
for this is that: i) t he Newtonian results show that the g-mo des of the ocean have frequencies 



that scale as ( Bildsten. Ushomirskv and C utler 



r-modes, would have frequencies of order Q ((Lockitch. Friedman and Anderssonl 12003). In 



1996 1), and ii) inertial modes, like the 



other words, at least for these classes of modes, our equations should be consistent. 



3 THE TRADITIONAL APPROXIMATION 

We now wi sh to analyse the various oscillation mod es of a neutron star ocean. To do this we 



follow e.g. 



Bildsten. Ushomirskv and Cutlerl (|1996[ ) and introduce the traditional approxi- 



mation. As in the Newtonian problem, this corresponds to 

• neglecting the Coriolis term in the r-momentum equation, 

• assuming that the mode is mainly horisontal. 

In our case, the first assumption corresponds to neglecting the H 2 term in equation (J26|) . 
while the second implies that we should neglect W in equation (|2"K|) . 

With these simplifications equations (J27)) and (J28|) can be solved for H\ and H 2 - This 
then gives us 

dpi 

—i{a + mVt) (p + p)r 2 sin 0(1 — q 2 cos 2 9) Hi = sin 9—— + mq cos 9pi (29) 

o9 
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and 



(a + m£l){p + p)r 2 sin 9(1 — q 2 cos 2 9)H 2 = 2q sin 9 cos 9— — + mp\ 

o9 



Here we have introduced 



Q 



2(n 



a + mQ 

By using the above results in (J2H)) we find 



(30) 



(31) 



i(cr + mtt)p l + (p+ p)e v W '- 



(P + p)[^ 



where we have defined the operator 

l-/i 2 d 



L = £- 



>^ 2 
v + - 

r 



m 



e u W 



le 



o- 



mVt)r 2 
mq(l + q 2 fi 2 ) 



L( P i), (32) 



(33) 



dfj, \1 — q 2 \i 2 dfx J (1 — /i 2 )(l — q 2 p 2 ) (1 — q 2 \i 2 ) 2 
with fi = cos 9. Notably, this has ex actly the same form as the angular ope rator in the 

Newtonian problem, cf. equation (7) in Bildst en. Ushomirskv and Cutler! (J1996). One has to 



be careful here, though, because in the relativistic case L depends on r through q (which in 
turn depends on ou(r)). However, as long as we are interested in oscillations in a thin radial 
shell we may take the frame-dragging to be a constant. 

Aiming to write the final equations in a form reminiscent of the Newtonian results, we 
relate the displacement component £ r to the perturbed velocity component via 

6u» = q» v £ut\ (34) 
where the projection q^ u was defined earlier, and where £ u is the Lie derivative along 

This means that 

e v W(r, #)e*(°" t+m< ^ 



i(a + mQ) 
e v H 1 (r,9)e i{r7t+m4 

i(a + mfl) 
e u H 2 (r,9)e i{at+m4 



(35) 
(36) 
(37) 



i(a + mVt) 

We have used the inherent gauge freedom to set the temporal displacement to zero, i.e. 
= 0. We also relate the Eulerian perturbations in density and pressure in the standard 



way; 
Ap 

r x p p + p 



A Pl _^ Pi _ Pi +£ A 7>;ai 



Tip 



P + P 



where is the adiabatic exponent Ti 



Pi 



dlnp 
<9 hip 



P + P 



Pi-C 



P + P , 

iip 



(38) 



evaluated at constant entropy. 



Given these various relations we can write equation 



-r- 



pi 



Tip \Tip 



P' 



X 
2 



as 
pe u 



(cr + mVt) 2 r 2 {p + p) 



(39) 
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Finally, the r-momentum equation (|26[) reads (after neglecting the Coriolis term) : 



(P + PY 



[a + mVtf 



-A 



e r =Pi + 2f^[p+(l + ri)p]pi, 



where we have defined the relativistic Schwarzschild discriminant as 

e -A/2 



P + P 



P 



P + P / 
Tip 



(40) 



(41) 



For later convenience, we also introduce the relativistic analogue of the Brunt- Vaisala fre- 
quency as 

N 2 = - V -A (42) 



4 THE NEWTONIAN LIMIT 



As a consistency check of our derivation we want to compare our final formulas (|39|) and (|4U|) 
to the corresponding results in the Newtonian case. To facilitate this comparison we need to 
transform our expressions, which are determined in a coordinate basis, into an orthonormal 
basis. T his basis should be the relativistic analog ue of the flat space-time orthonormal basis 



used by 



Bildsten. Ushomirskv and Cutlerl (|19961 ): 



(43) 



where the subscript N denotes a Newtonian vector. 

In the corresponding relativistic basis, the spatial components of the displacement 
vector are determined from 



C = <\ 
^ = rsin0£« 



(44) 
(45) 
(46) 



In this frame we can identify the components of the displacement with their Newtonian coun- 
terparts. Moreover, we can ar gue in favour of the validity of th e trad itional approximation 



in a way analogous to that of 



Bil dsten. Ushomirskv and Cutlerl (199 



In the orthonormal frame equation (J2T?)) becomes 



-I/+A/2 



{a + mfi) 2 + A 



C = Pi + 



V 



P + P 



+ 1 Pi 



+ie~ u {a + mQ) [(Q - u) {rv' + 2) + rJ\ sin B£*. 



(47) 
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We can safely neglect the Coriolis term (the last term of the right-hand side) as long as 



N 2 Q > e~ v (a + mVL) - u) {rv' + 2) + rJ\ sin0f£| 



(48) 



or equivalent ly 

|iV 2 | > (a + mn) (rv' + 2)+ru']\ 



So 



So 



(49) 



Since we are assuming that the mode is horisontal, it is far from obvious that this inequality 
will hold. However, for a neutron star ocean the order of magnitude of ruj' can be estimated 
from the value at the surface. From equation (jl4j) . we see that 

6IQ 12MQ 



ru' « Ru'(R) = 3u(R) 



R 3 



5R 



(50) 



12 



and since M/R « 1/5 for a typical neutron star it follows that ru' ~ 25^' ■'■ n a siixiilai" way, 



2e A M 2M 



r 2 R 2 

which means that rv' « 2/5. Combining the various terms the inequality becomes 



(51) 



\N 2 \ > 1.5aQ 



(52) 



According to iBildsten and Cutler! (|l995T ). the thermal buoyancy in the deep ocean leads 
to iV ~ 1—5 kHz. Using the slow- rotation solutions one can also argue that |£o/£ol ~ 
10 4 m/10 2 m = 10 2 . This means that, for mode frequencies below (say) 10 Hz our approxima- 
tion remains valid for rotation frequencies up to a few hundred Hz. Given that the highest 
observed spin rates in Low-Mass X-ray Binaries are about 600 Hz, our formulas are likely 
to be useful for studies of low frequency ocean waves in most accreting neutron stars. 
In the orthonormal frame our two final equations are 



-u+X/2 



V l e v-\I2 ■ 

(cr + mQ) 2 H A 



So 



v P+P 



1 pi- 



(53) 



and 



pT 1 r ' q ° 



e x/2 p 1 



e u+x/2 p 



-L 



(54) 



pTi (p + p)(a + mVt) 2 r 2 

We want to verify that these equations limit to the expected Newtonian ones. To check this 
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we note that the Newtonian limit corresponds to 

|$| <C 1 where $ = —M/r is the Newtonian gravitational potential , 
p <^ p 
u -> 



v 



2$' 



1 - 

5 

2M 



2M 



1 + 2$ « 1 



2<? 



where g is the gravitational acceleration. Using these simplifications, letting — > 

as well as introducing the scale - height /i = p/pp, equation (J53|) yields equation (5) of 



Rildsteu. Ushomirskv and Cutler! ()199ffl: 



p(a 2 r - N 2 )Ck = A + 



V 



(55) 



where a r = a + mQ is the mode frequency in the rotating frame. At the same time, equation 
reads 



r x p 



+ 



PS 2 N 



(56) 



which is equation (6) of Bildst en. Ushomirskv and Cutler! (J1996J) (where the term 2^ r N / r was 



omitted due to a misprint). 



5 THE OCEAN G-MODES 

We now want to assess the magnitude of the relativistic effects on the ocean g-modes of a 
rotating neutron star. In principle, we could do this by solving ()53|) and (|54p numerically, 
thus finding the modes of oscillation for a given model of the ocean. It may, however, be more 
instructive to extract the required information through an approximate analytic calculation. 
As this is the route we have chosen, we first simplify our equations somewhat by assuming 
that 

P « P, (57) 

e" A = e u « 1 , 58 

r 

v' « __ = 2e A / 2 ^, (59) 
2 J 
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1 r 



0.1 1 10 

2(Q.-(o) 

q = — 

(a+mQ.) 

Figure 1. An illustration of the relation between the numerically determined angular eigenvalues a and the spin parameter q. 
The eigenvalues for the g-modes asymptote to ctRJ (2n — l) 2 q 2 for n = 1, 2, 3... while the trapped Kelvin modes correspond to 
a fa m 2 for large q. 



These approximations should be quite accurate in the low-density ocean of a neutron star 
We also denote the "angular eigenvalues" by a, i.e. we assume that 

( Pi 

p j V p 

In Figure H we show a sample of numerically determined ang ular eigenvalues a. These re- 



LI* 



(61) 



suits w ere calculated using a method similar to that employed by 



Bildsten. Ushomirskv and Cutler 



(|l996l ). These results clearly illustrate the large-g asympt otic behaviour. It should be noted 



that t h e eigenvalues can be classified in different ways, cf. iBildsten. Ushomirskv and Cutler 
(119961 ): iLee and Said ()1997l ). For example, in the q — > limit we can associate each solution 
with the I and m of the single spherical harmonic eigenfunction. Note that our data corre- 
sponds to I = 1 and I = 2 and the 21 + 1 solutions that correspond to the various permitted 
values of m are easily distinguished in the figure. Alternatively we could label the modes by 
their symmetry with respect to the equator (it depends on whether / + m is odd or even). 
It should be noted that, for large q, the first three g-mode asymptotes then correspond to 
odd, even and odd eigenfunctions, respectively. 



Having defined the angular eigenvalue a, we obtain the two equations 



-e-< 



2M 



-1/2 



Pi 
Tip 



2M 



-1/2 



P9 
Tip 



2M 



-1/2 



api 



p(a + mVt) 2 r 2 
(62) 
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and 



9 



1 



2M 



-3/2 



{a + mQ) 2 + 1 - 



2M 



-1/2 



gA 



Co=P\ 1 



2.V\ 1/2 p# 



Tip 



Pi- (63) 



These equations simplify considerably if we introduce new variables 77 = r 2 i/j^ r and 
= pi/if), where 



ip = exp 

and the "sound speed" follows from 



In terms of these new variables we get 



2.U\~ 1/2 g 



dr 



-7] 



and 



1 - 



1 



2M 



-1/2 



1 



Ftp 



2M 



a 



2M 



-1/2 



(a + mVtf 



r J (a + mVt) 2 r 2 p 
2M\ 



,2 1 /,2 



1 



- gA 



r I 



9V 

2,/,2 



r 2^/, 



(64) 



(65) 



(66) 



(67) 



If we now assume that the radial dependence of the perturbations is described by rj ~ 
^ exp(zfcr), where k is the radial wavenumber, we readily arrive at the dispersion relation 



2M 



1 )gA 

r 



9 



2M 



a 



r J (a + m£l) 2 r 2 



(68) 

We now use the Brunt- Vaisala frequency iV 2 = —gA and an analogue to the standard Lamb 
frequency; 



-2 = ac s 



Implementing a dimensionless corotating oscillation frequency 

a + mQ 



K 



n 



we get the final result 



k 



O 2 
c 2 /t 2 



2M 



-1 



K 



2M\ ( N 



K 



2M\ ( L 



r J \Q 

In o rder to estimate t he g-mode frequencies we recall that (cf. 
(1996); Townsendl (2003)) 



n 



(69) 



(70) 



(71) 



Bildsten. Ushomirskv and Cutler 



a w (2n - l) 2 g 2 = (2n - 1) 



2(1- w/fi) 



n 2 



for large q fsee lBildsten. Ushomirskv and Cutlerl (|1996[ ) for a discussion of the interpretation 



n = 1, 2, 3, ... 



(72) 



of the "quantum number" n). This behaviour is apparent in Figure Expanding equation 
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(J72*Jl for low frequencies (that is, for k « N/Q and k « L a /Q) we find 



^^-^V-T^-'-t) {m) • (73) 

Even though this does not actually determine the mode-frequencies it allows us to estimate 
the magnitude of the relativistic effects. We see that there are two main effects: The obvious 
one is the gravitational redshift factor, which for a typical neutron star with M/R « 1/5 
decreases the frequency by about 10% compared to the Newtonian case; 



1/2 / 2M\ 1/4 / iV N 1/2 



2M\ l/A M 



— « 0.9. (74) 



i? / 2i? 

The rotational frame-dragging adds a further 10% to the decrease since 



2J 2IQ AMQ 

to 



1-^ 



1/2 

« 0.9. (75) 



R 3 R 3 5R 

Thus, we would expect the relativistic ocean g-modes to have frequencies roughly 20% lower 
than their Newtonian counterparts. 

It is easy to argue that these relativistic corrections must be included if one wants to draw 
the correct conclusions from observed oscillations. For example, comparing our relativistic 
result to the corresponding Newtonian one (where u and M/R are both taken to be zero 
in (J73J)) we see that already for n > 2 would we get kgr( u + 1) < K Ncwton(^)- This means 
that the Newtonian formula would predict that the wrong mode was observed. This could 
be crucial if the aim is to study the inverse problem to put constraints on the physics of 
neutron star oceans. 

As is clear from Figure ^ there is a distinct class of angular eigenvalues which do not 
exhi bit the sca l ing d iscussed above. These solutions correspond to trapped Kelvin waves, 
and Townsend ( 20031 ) has shown that for large q one gets 



2 2mq 2 

a = m ~ m (76) 

2mq + 1 K J 



Following the same procedure as for the g-modes, we now find that 

2M\ 1/2 N 



K 171 



(77) 

R J ttkr K J 

The derivation of this relation is valid as long as the radial wavenumber k is taken to be large, 

i.e. for short-wavelength waves. Otherwi se we would not have k « N/Q and k « L a /Q. 



Now comparing to the Newtonian result ( Townsendl 
the frequency by about 20%, since 



2003) one finds that the redshift reduces 



2M\ 1/2 M , . 

1-— J wl--«0.8. (78) 



compared to the Newtonian case. 
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6 OCEAN INERTIAL MODES: THE R-MODES 



i whether 



Townseri' 



Given the relativistic version of the equations of motion it is interesting to as 
we ca n estimate the frequencies of r-modes in the ocean. As discussed by eg. 
( 2003 ) the r-modes are not well represented by the traditional app roximation. However, 
it is straightforward to follow iProvost . Bcrthomieu and Rocca (1981) and assume that the 
following scalings apply (see also 



Lockitch. Friedman and Andersson 



(2QQ2()): 



Pi ~ pi ~ Q 2 , 

The relativistic continuity equation then immediately reduces to 



im sin 8H2 = 0, 



(79) 



and a second equation relating Hi and H 2 can be obtained by combining the 9- and <p- 
momentum equations. This leads to 



2i(a + 2mQ — muj) sin 9 cos 9H 2 + i(a + mQ) sin 2 6 



dH 2 
09 



+ 



[m{a + mO) + 2(0, - u)(cos 2 9 - sin 2 0)]H X + 2(fi — uj) sin 9 cos 9 



dHi 
~~dT 



0. 



(80) 



After some straightforward algebraic manipulations one can show that these two equa- 
tions will be satisfied if 

1 



#1 



sin 9 



and 



a = — m£l + 



Pi m (cos9) 



2m(Q — u) 



(81) 



(82) 



1(1 + 1) ' 

It should, of course, be noted that this derivation only works as long we can take the frame- 
dragging, u, to be a constant, i.e. for oscillations in a thin shell. Anyway, we conclude that the 
frequencies of the ocean r-modes are sensitive to the surface value of the frame-dragging. The 
size of the effect can be appreciated if we take the ratio between the relativistic co-rotating 
frequency kqr and its Newtonian counterpart ^Newton- This immediately leads to 

kgr 1^1 84 (83) 

^Newton ^ 5i? 

Thus, the relativistic frame-dragging lowers the r-mode f requencies (in the rotating frame) by 



about 15%. This is in accordance with the resu 



(|2003|) and 



Abramowicz. Rezzolla and Yoshida 



t s of e .g. 
(20001). 



Lockitch. Friedman and Andersson 
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Finally, it is worth pointing out that the calculation of general inertial modes is much 
more involved since one must account fo r the coupling of many spherical harmonics, see 



Lockitch. Friedman and Anderssonl (j2003l ) for a recent discussion. In particular, such a cal- 



culation may be needed for the fastest spinning neutron stars since the g-modes will also be 
dominated by the Coriolis force when Q » N. 



7 CONCLUDING REMARKS 

We have shown how the so-called traditional approximation can be generalised to general 
relativity. This leads to a set of equations that ought to prove valuable for future studies 
of waves in the oceans of rotating neutron stars. The derivation of these equations required 
several assumptions in addition to those made in the Newtonian case. First of all, we have 
to make use of the Cowling approximation. That is, we neglect all metric perturbations. 
This may seem a drastic simplification, but it is a reasonable approximation for waves in 
the relatively low-density neutron star ocean. After all, such waves are not expected to lead 
to appreciable gravitational-wave emission. Secondly, we have to assume that the frame- 
dragging Co is constant, which restricts us to consider thin shells. Again, this should be a 
good approximation for neutron star oceans. Finally, and perhaps most interestingly, we 
have to neglect some terms of post-Newtonian order from the perturbed relativistic Euler 
equations, an approximation which is valid for low-frequency oscillations. Having made these 
various assumptions we can separate the eigenf unctions into a radial and an angular part. 
The angular eigenvalue problem is formally identical to the Newtonian one: it corresponds to 
solving Laplace's tidal equation, which now depends also on the dragging of inertial frames. 
It should be emphasised that, since the two problems are identical we can use the results 
for the Newtonian angular problem also in the general relativistic case. 

Our calculation prepares the ground for detailed calculations based on sophisticated 
models of the neutron star ocean. Given the eigenvalues to Laplace's tidal equation, we 
need to solve the two radial differential equations (jB^j) and ()54|) . together with the relevant 
boundary conditions. Once one provides a description of the physics of the ocean, the solution 
of these equations should be straightforward. Here we opted to discuss the effect that general 
relativity has on various classes of ocean modes in a less quantitative way. We simplified 
our equations to the conditions that should prevail near the surface of a neutron star. This 
allowed us to derive a dispersion relation which we used to deduce to what extent the various 
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classes of modes are affected by the gravitational redshift and the rotational frame-dragging. 
Thus we have shown that the mode frequencies are typically lowered by roughly 20%. This 
is not particularly surprising, but it should be noted that the overall effect is due to different 
combinations of the frame-dragging and the redshift factors for the different kinds of modes. 
We also discussed (briefly) the effect that general relativity has on the ocean r-modes. 

Our estimates provide a strong argument in favour of using the derived formulae in more 
detailed (numerical) studies. In particular, we believe it is clear that one must account for 
the relativistic effects if the aim is to identify the individual ocean modes present in observed 
data. 
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